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(54) Constrained shortest path routing method 

(57) A real-time method for routing subject to an ac- 
ceptable delay constraint between nodes in high-speed 
data networks, such as PNNI protocol networks, uses 
an illustrative two-phase algorithm employing Dijkstra's 
algorithm at each phase. In an illustrative first phase, 
the Dijkstra SPF algorithm is used in seeking the short- 
est cumulative delay from the destination to the source, 



thereby generating cumulative delay labels from a node 
/to the destination node k. The delay results are then 
employed in the second phase, where the Dijkstra SPF 
algorithm is illustratively employed for determining ad- 
ministrative weight (AW) as the link metric subject to 
modification in accordance with results obtained in the 
first phase. 
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Description 

Field of the Invention 

s [0001] The present invention relates to methods and systems for routing data paths in high-speed data networks. 
More particularly, the present invention relates to methods and systems for selecting preferred routing paths in such 
networks employing private network node interface (PNNI) protocols. 

Background of the Invention 

10 

[0002] High capacity data networks, such as those using high-speed asynchronous transfer mode (ATM) switches 
offer promise not only of high-speed data delivery, but also a variety of quality of service (QoS) guarantees. An important 
element of establishing paths interconnecting nodes in a such high-speed data networks is an efficient, reliable routing 
mechanism. 

is [0003] Recently, the ATM Forum has standardized on certain PNNI protocols. As defined, these protocols provide 
for signaling and routing protocols that permit connection setup and release with multiple QoS parameters between 
nodes. These protocols also provide for exchange of topology state information between nodes. See further, "Traffic 
Management Specification," ATM Forum, 95-0013, R10, Feb. 1996; and Dykeman, D., and M. Goguen, "PNNI Spec- 
ification," Ver. 1 .0, ATM Forum, AF-PNNI-0055.000, May, 1996. 

20 [0004] Many implementations of the PNNI protocols have been proposed, with a variety of performance and QoS 
features. An important consideration in many networks, such as those using the PNNI protocols, is interconnecting 
desired nodes while employing a minimum of network resources. In particular, it is desired to interconnect nodes using 
the shortest interconnecting paths. Additionally, it has become ever more important - especially for many time-depend- 
ent applications (such as multimedia applications) - to seek to achieve low or minimum delay time for transmission 

25 between selected nodes. QoS criteria based on such reduced or minimum delay are therefore of great importance for 
critical applications, and provide important competitive differentiation for service providers. 

[0005] There are many solutions for selecting the shortest path (or minimum defined cost) between selected nodes, 
even in a network having a large number of interconnected nodes. Prominent among these shortest path algorithms 
is the well-known Dijkstra algorithm. Solutions based on the Dijkstra and other algorithms may often be modified to 

30 permit the inclusion of certain conditions relating to various QoS criteria. Experience with Dijkstra algorithm solutions 
and corresponding coding implementations have proven reliable in many cases and have been adopted for realtime 
network implementations. It is desirable, therefore, to reuse existing code and maintain functional compatibility between 
delay-constrained solutions and existing shortest path implementations. As is well known in the art, however, network 
solutions of the shortest path problem subject to an additive delay condition is a so-called "NP hard" problem that 

35 proves intractable in real-time network contexts using current implementation technologies. 

Summary of the Invention 

[0006] The limitations of the prior art are overcome and a technical advance is made in the network routing arts in 
40 accordance with illustrative embodiments of the present invention, as described in the following detailed description. 
[0007] In accordance with an illustrative embodiment ol the present invention, determinations are made at a source 
node in a network of the short weight paths to each other node in the network, subject to satisfying an acceptable delay 
constraint. Advantageously, these determinations are performed in a two-phase method employinga modified Dijkstra's 
algorithm at each phase. 

45 [0008] In an illustrative first phase, the Dijkstra SPF algorithm is used in seeking the shortest cumulative delay from 
the destination to the source, thereby generating cumulative delay labels from a node; to the destination node k. The 
delay results are then employed in the second phase, where the Dijkstra SPF algorithm is illustratively employed for 
determining administrative weight (AW) as the link melric subject to modification in accordance with another aspect of 
the present invention. More specifically, in an illustrative embodiment it proves advantageous to (i) label a node / with 

50 cumulative AW from the source node, (ii) track cumulative delay (D s j) from the source during forward AW labeling, and 
(iii) labeling neighboring nodes / of a permanently labeled node / only if (D s { + d j k + Sjj) is less than the specified end- 
to-end delay constraint (where Sj j is the delay for the link between nodes / and /). 

Brief Description of the Drawing 

55 

[0009] The present invention as described in illustrative embodiments in the following detailed description and in the 
attached drawing, wherein: 

[0010] FIG. 1 shows an example network for illustrating an illustrative embodiment of the present invention. 
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[0011] FIG. 2 is useful In understanding Phase 1 of an Illustrative embodiment of the present invention as applied to 
the example network of FIG. 1. 

[0012] FIG. 3 is useful in understanding Phase 2 of an illustrative embodiment of the present invention as applied 
to the example network of FIG. 1 . 

Detailed Description 

[0013] FIG. 1 shows an illustrative network of nodes corresponding, for example, to the nodes in an ATM network. 
Each source (originating) node includes a control processor for receiving network node state information (as is well- 
known) and determining a low-cost or near-optimum shortest path analysis, subject to a maximum cumulative delay 
constraint. Costs are illustratively measured, in part, in terms of administrative weight (AW). Additionally, in an accord- 
ance with another aspect of the present invention, costs are also measured, in parti in terms of link delay. Initially, the 
present discussion proceeds in terms of a single source seeking a low-cost routing to a single destination {unicast 
operation). Single-source, multi-destination (multicast) operation will be discussed below. 

[0014] Routing is performed, for a given caH, at the source or on behalf of the source by an associated separate or 
shared processor. In accordance with PNNI protocols, link state parameters communicated between the nodes are of 
two types: link metrics and link attributes. A link metric is a state parameter for a link, which parameter is a combination 
of all link parameters along a given path. Thus, such a link metric can indicate whether a particular path is acceptable 
for a given connection. A link attribute, on theother hand, is a link state parameter that reflects the state of a given link 
in a path. Path delay is a link metric suitable for use in path weight determinations in accordance with one aspect of 
the present invention (rather than AW). Path delay is taken as the sum of the delay of all of the links in a path. Received 
path delay information from other network nodes is advantageously stored in table form at a node processor. 

Model 

[0016] It proves useful to more precisely characterize network routing solutions in terms of a network model. For this 
purpose, let G = (V t E) represent a directed graph representation of a network. We assume that there are n nodes and 
m arcs in the tree. Let s e V represent the source and f represent the destination. Let (/, /)££ denote a directed edge 
(arc) in the graph. Let 0 represent the weight of arc (i t j) and dg > 0 represent the delay incurred along-arc (/; j). 
Let P(i, k) represent the set of simple paths from node / to node k. The weight of a path PC P(/,/r), denoted by w(P), 
is defined as the sum of the weight of the arcs in the path. The delay of a path P G P(l A), denoted by d(P), is the sum 
of the delay of the arcs in the path. The objective of the routing problem is to find the path of minimum weight from 
"source s to destination t such that the delay of that path is less than some pre-specified value D, i.e., 



in w(P) 



subject to 

d(P) <;D. 

As noted above, this problem is NP-hard. Next, we develop a heuristic algorithm that seeks to minimize the weight of 
the path while avoiding violating the delay constraints. We show that the algorithm always terminates with a feasible 
path (with respect to the delay constraints) if one exists. We will also show in the sequel that implementations of the 
present inventive algorithm obtain paths of low weight. 

Algorithm Description 

[0016] In this section, we describe the algorithm that is run at each source. We first describe the algorithm for a 
particular destination. This algorithm is then used as a subroutine to compute short weight path to all the destinations 
from a given source. 

Single Sink 

[0017] Let s represent the source node. In this section, we determine a short weight path from s to a given destination 
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node r. The processing of the algorithm at a node processor advantageously proceeds in two phases. Phase 1 of the 
algorithm involves running the well-known Dijkstra algorithm for paths from the destination to the source using the 
delay labels as the length of the arc. The distance labels (in delay) that result from the execution of Phase V are used 
in Phase 2 of the algorithm. We now give the algorithm for Phase 1 . 
S PHASE 1 (s, t, d) 

• INITIALIZATION 

- U= V, W = <J>, 6,= 0, 5 y =ooV y .* t 

10 

• ITERATIVE STEP 

- k = Argminjet/8y 

is - U= U\[k\ and W^Wkj {k}. 

- For each /€= U, {j, k) G E, if 6*+ d jk < Sythen 5 y = 6* + d jk 

• TERMINATION TEST 

20 

Stop if s e W. 

[0018] Note that bj represents the shortest path length from node y to the destination t. If 8 S > Dthen there is no 
feasible solution to the problem and the algorithm terminates. If 8 6 <> 0, then Phase 2 of the algorithm, which is a 
25 modified version of Dijkstra's algorithm, is executed from the source s. 

[0019] For the ease of the algorithm description, we define the generalized weight of an edge as 



30 A (Wfrdi,) (1) 



The generalized weight of a simple path PG F\i, k) is defined as v (P) • (w (P), d (P)). 

[0020] In Phase 2 the predecessor of each node is maintained in n [i], V/G V, and atwo-tuple label (go, is maintained 
35 for each node, where eo is the accumulated weight from the source s along a feasible path and y is the delay from the 
source s along the same path. A lexicographic order is defined for these two-tuples, namely, 
Definition 1 Define 

if 

• ©/ < toy, or 

45 • (Oj^&jandy^yj 

[0021] Phase 2 of the algorithm is: 
PHASE 2 (s, f, w, cO 

so . INITIALIZATION 

1: U= V, W - 0 t co 3 = \\f s =0, <oj = \\fj = °°Vj* s, 

• ITERATIVE STEP 

ss 2: k = Arg min yeL , (a>y> Vy) 

3: U=U{K}and lV=WWty 
4: 4.1: For each y S U, (k, j) G E 
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4.2: if (<o k + y k + d^) < (a^ y y ) 

4.3: MV/c^w +*/*£>) 
4.4: coy = <o k + 
4.5: v y = \(f| C + c(i ? 
5 4.6: jtfl] = k 

• TERMINATION TEST 
5: Stop re IV. 

10 [0022] Note that if line 4.3 is taken out, the above algorithm is just Dijkstra's algorithm with respect to the generalized 
weight v. In other words, it is Dijkstra's algorithm with respect to the weight with the delay label \j/ used as the tie- 
breaker. Since all the arguments used in proving Dijkstra's algorithm still apply when we use the generalized weighty 
defined in 1 , we have the following lemma. 

Lemma 1 1f we run Dijkstra's algorithm on G with respect to the generalized weight v, then at termination, 

is 

20 [0023] The following theorem summarizes the properties of our algorithm. 

Theorem 2 Phase 2 terminates with a path that satisfies the delay constraint. In addition, if there are shortest paths 
with respect to weight that satisfy the delay constraint, the algorithm will find one. 

Proof. 

25 

[0024] First we show that when 6 S £ D, Phase 2 will find a path that satisfies the delay constraint. 
[0025] From line 4.3 and 4.4, one can see that the weight label of a node /, a>j, assumes a finite value only when 
node / is on a feasible path (/.e., there is at least one path from /to t that satisfies the delay constraint). In addition, line 
3 removes one element from set U \n each iteration. Therefore, if the algorithm fails to find a feasible path, we must 
30 have ©,= Vf e Uat some point during the execution of the algorithm. This is not possible due to the following reason: 
Whenever a node is included in the set W there exists a path from the source to the destination through that node that 
satisfies the delay constraint. Since s E Wand t G U (otherwise the algorithm has succeeded), 3k e WJG U, s.t. 
edge (k,j) is on a feasible path. When k was included in WVr\e value of Wj would have been updated to a finite value 
and hence a contradiction. 

35 [0026] To prove the second statement of the theorem, first consider running Dijkstra's algorithm with respect to the 
generalized weight v. When it terminates with a simple path from s to t, denoted by P*(s,0 =(s,/ 1t .... ^ = 0. each node 

on this path has a label (<o, v k = 1 ko. We know from Lemma 1 that P*(s, t) is the shortest path with respect 

to weight among all the paths from s to f, and if there are more than one such shortest paths, P*(s,f) has the smallest 
delay among them. Therefore, if there are shortest paths with respect to weight that satisfy the delay constraint, P*(s, 

40 t) is one among them. 

[0027] Let's assume for now that no other shortest path with respect to weight has the same delay as P*(s, f). We 
then use mathematical induction to prove th^t upon termination of Phase 2, a) the label for node i k on P*(s, 0 » s still 
(co //c , y ik )' t and b) the path from s to t will be found to be P* (s, /)• 

[0028] Since P* (s, 0 is a feasible path, the condition on line 4.3 will be satisfied when are (s, /,) is explored. By 
45 Lemma 1, arc are (s, a,) has the smallest generalized weight v rf|l among all the paths from sto /,, so node ^ will be 
permanently labeled as (<o/,,y/i) =(w sir d sil ) and its predecessor will be s. 

[0029] Note that it is possible that there are other paths from s to J, whose generalized weight is the same as (w siv 
d si) ), in which case node t\ may be permanently labeled not from s but from the last node on one of these paths, say 
the last node of P'(s, But since the value of the label for /, will be the same, and s will still be on the predecessor 
50 tree of / 1f possibly with some intermediate predecessor nodes, finding path P\s, i^) is no different from finding (s, /,) 
as far as weight and delay are concerned. Therefore, we will ignore such cases without loss of generality in the rest 
of the proof. 

[0030] Since node /, is directly connected to s, it acquires its permanent label when s is taken out of set U, which 
happens during the first iteration of the algorithm. Because the label values do not increase and (o> //f , y ik ) £ (co /fcf v V/> +1 ) 
ss for k= 7,... , Aq - 1 , we know that if (<o iv y h ) < (o fel \{/ fe ) (or equivalently (w hi2 , d hi2 ) > (0,0)) then node /, will be removed 
from set U before i 2 , i^. If {<o h , y h ) = (o^, y fe ) (i.e. (iv /lfel d hiz ) = (0,0)), then it is possible that ^ is removed before 
/, from set U. But this simply means that we have found another path from s to ^ with the same weight and delay as 
P*(S,fe) - (&, /,,£)- Again, we dont distinguish such cases from that of t, being removed frorn set U before fe, for-these 
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paths are equivalent as far as weight and delay are concerned. For the rest of the proof, we assume (o) /v <(<*> //f+1 , 
Vfc^) without loss of generality. 

[0031] Assume for the purpose of induction that the nodes f, i k have been moved to set Wand permanently 

labeled as (co /v y,,),..., (co ik , Further assume that the path from s to i k \$ found to be (s, i k ) and is removed 

s from set Ubefore nodes i kn , ...J^. We now show that node willbe the next on the path and will be labeled (<o /it+1 , 

[0032] When node i k is removed from set U, arc ) is explored. Since the path from s to i k is part of P*(s, f) and 
i M is on P*(s, 0, condition on line 4.3 is satisfied. Hence i k +, will be labeled as (<o ik , y ik ) + (^y^d^) = (co /Apf1 , (¥/* + ,), 
and its predecessor will be i h Because (w^^y^,) is the smallest label by Lemma 1, node k +1 will not be relabeled. 
10 in addition, i k ^ will be removed from set U before nodes i^ +2 , .... i^ because it acquires its permanent label when node 
i k is removed from set U and this label is smaller than those of i k +&...j k(r This concludes our induction. 
[0033] Finally, if there are other shortest paths with respect to weight that has the same delay as P* (s, f), then it is 
possible that such an alternative path is found instead of P*(s, Q. But these paths are equivalent as far as weight and 
delay are concerned. 

15 

Multiple Sink Case 

[0034] The objective in PNNI is to determine the shortest path from source s to all other nodes. This is done by 
executing Phase 1 and Phase 2 from the source to each sink. However, when the algorithm is executed to determine 
20 a short weight path from the source to the destination tj, it might give information about short paths from the source 
to other destinations. In particular, for all destinations in the set F, the optimal solution to the constrained shortest path 
problem is known. This is exploited in the multiple sink case. The algorithm multiple sink is given below: 

MULTIPLE DESTINATIONS(s) 

[0035] 

• INITIALIZATION 

30 - G=VF=0M=~. 

• ITERATIVE STEP 

- Pick re G 

35 

- PHASE 1 (s, t, d) 

- PHASE 2 (s, f, w) 
40 . TERMINATION TEST 

- Stop if G G 0 else go to the ITERATIVE STEP 

[0036] Returning now to the example network of FIG. 1 , we consider the application of phase 1 processing described 
45 above. As illustrated in FIG. 2, phase 1 comprises determining the smallest delay labels (dj k ) for the delay from a node 
/to the destination /fusing the reverse Dijkstra SPF algorithm. It is assumed that the Source-to-Destination delay is 
constrained to be less than or equal to 10 units. The link parameters are presented as (link AW, link delay). So, for 
example, for the topmost path in FIG. 2, the link (AW, delay) tuple is (1 ,4). 

[0037] In FIG. 3, the results of phase 1 processing are applied in phase 2 processing. In particular, we again run the 
so Dijkstra SPF, this time from the source. AW is used as the link metric for the "shortest" path determination, but (i) each 
node is labeled with cumulative AW sum from the source, (ii) cumulative delay from the source (D s{ ) labels of a neigh- 
boring node j of a permanently labeled node /are updated only if (D s i + dj k + 5j tj ) meets the end-to-end delay constraint 
(10, in this example). It becomes clear that the path along the upper periphery of the network representation would 
have a delay of 11, which exceeds the delay constraint. The path from the source through nodes 6 and 3, however, 
ss has a total path delay of 9. The tuple inside the node circles reflects the sums of link AW and the sums of link delay, 
respectively. 

[0038] From the foregoing, it will be appreciated by those skilled in the art that currently available Dijkstra algorithm 
technology, including existing program code, can be adapted in accordance with the present inventive teachings to 
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allow real-time routing at nodes of a data network in response to link data provided in the normal use of PNNI protocols. 
[0039] In use of the present inventive methods and systems employing disclosed algorithms and methods, node 
processors or alternative processor facilities will find one or more feasible paths through a network that satisfy end-to- 
end delay constraints - if at least one such path exists. 

s [0040] If the AW shortest-path automatically satisfies the delay constraint, that path will advantageously be chosen 
using the present inventive teachings. While the above-described methods indicate the use of the Dijkstra algorithm 
n + 1 times for the one-to-all shortest path solution, those skilled in the art will recognize that in appropriate circum- 
stances the computational load can be reduced by exchanging delay-to-destination labels dj k among nodes. While 
the present detailed description has proceeded in the context of PNNI protocols, the present invention will prove useful 

io in a variety of network types using a variety of network protocols. 

[0041] Other modifications and extensions, all within the spirit and scope of the attached claims, will occur to those 
skilled in the art in light of the present inventive teachings. 



is Claims 

1. A method for choosing a path in a network from a source to a destination, the network comprising a plurality of 
nodes, each node connected to a plurality of other nodes by links, each link having a weights reflecting at least 
first and second link parameters, paths connecting said source and said destination comprising a plurality of links, 

20 said paths being subject to a constraint on the cumulative value of said second parameter for links in a path, the 

method comprising the steps of: 

a. for paths from said destination to said source, performing a minimum-weight path identification based on 
said second parameter for paths in said network, thereby to generate labels identifying weights for said second 

2S parameter for each node /to destination k t 

b. for paths from said source to said destination performing a minimum-weight path identification from said 
source using said first parameter as a link metric, and 

30 c. eliminating any nodes that would cause the constraint on said cumulative value for said second parameter 

to be violated. 

2. The method of claim 1 wherein operation of said network is in accordance with PNNI protocols. 
35 3. The method of claim 1 wherein said first parameter is administrative weight (AW). 

4. The method of claim 2 wherein said second parameter is link delay. 

5. The method of claim 1 wherein said minimum-weight path identification for paths from said destination is aceom- 
40 plished using a Dijkstra SPF algorithm. 

.6. The method of claim 5 wherein said minimum-weight path identification for paths from said source is accomplished 
using a Dijkstra SPF algorithm. 

is 7. The method of claim 6 wherein the method further comprises repeating steps a, b and c for paths from said source 
node to each possible destination node in said network. 

8. The method of claim 7 wherein for destination nodes after a first destination node, path information from said 
source to one or more other destinations is used in determining acceptable paths to each said nodes after said 

50 first destination. 

9. The method of claim 1 comprising the further step of eliminating any paths in said mini mum -weight path identifi- 
cation for paths from said source for which the cumulative value for said first parameter exceeds a predetermined 
maximum. 

55 

10. The method of claim 1 wherein said step of performing a minimum-weight path identification from said source 
comprises (i) labeling a node / with cumulative AW from the source node, and (ii) tracking cumulative delay (D s j ) 
from the source based on values for said second parameter during said minimum-weight path identification from 
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said destination. 

11. The method of claim 10 wherein said step of eliminating comprises labeling neighboring nodes j of a permanently 
labeled node / only if (D eJ + d iM + Sg) is less than the specified end-to-end constraint on said second parameter, 

s where 5g is the value of said second parameter for the link between nodes / and / and dj k is the value for said 

second parameter for the link between nodes yand ft 

12. The method of claim 10 wherein said second parameter is link delay. 
10 13. The method of claim 11 wherein said second parameter is link delay. 
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